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Abstract 

Electron transport is considered in high density fully 
ionized liquid metals. Ionic structure is described in terms 
of hard- sphere correlation functions and the scattering is 
determined from self-consis tently screened point ions. 
Applications to the physical properties of the deep interior 
of Jupiter are briefly considered. 



I. Introduction 


We are concerned here with the problem of calculating the 

resistivity of dense conducting fluids consisting solely of 

massive point ions and a neutralizing gas of interacting 

electrons. Several systems of physical and astrophysical 

interest are included in a calculation assuming the following: 

(i) that the density of the system is such that the electrons 

can be treated non-relativistically. If n £ is the electron 

-2 

density, this restriction can be stated as r » 10 where 

s 

r g is the usual linear measure of electron density 



r 3 a 3 
s o 



(ii) that the electron gas is degenerate. This is an implied 
restriction on the temperature, namely 


T « 


6x10' 


K 


s 

(iii) that the first Born approximation is adequate for the 

calculation of electron scattering cross sections from the 

ionic system. This condition is satisfied for r < \ 

(where +Ze is the charge on the point ion) and is discussed 

in detail in Appendix A. At lower densities (larger r ) 

s 

the validity of the results must be viewed with the caution 
normally attributed to low order calculations in liquid metals; 
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(iv) that the density-density correlation function (static 
structure factor) of the ionic- system can be approximated 
reasonably well by regarding the ions as an assembly of 
impenetrable spheres. In the presence of an electron gas 
(and with due account for the effects of exchange, correlation 
and the adiabatic response to ionic motion) the effective 
ion- ion interaction is characterized at short range by a 
steeply repulsive region, and at long range by a weak oscil- 
latory tail.^ At sufficiently high density (r « 1) 

the interaction between ions is expected to depart from the 
hard-core model and approach the simple screened interaction 
following from Thomas-Fermi theory (as used by Hubbard and 
Lampe 2 ); 

(v) that the contribution to the resistivity from electron- 
electron collisions can be neglected. So long as the electron 
system is highly degenerate this assumption is reasonable. 

In the following section we outline the basis of the 
calculations for the conductivity, and in subsequent sections 
estimate the melting temperatures of these fully ionized 
systems. The extensions to alloys are also discussed, and 
insofar as they apply the results are considered in the context 
of the physical properties of the deep interior of Jupiter. 



II. Calculation 


Within the adiabatic approximation we may write the 
resistivity of the dense ionized fluid of N ions in volume 0, 
as 


2 

P = m/n e t 
e 

where the transport relaxation time t is given by 


(1) 


1 

T 


« 


• 2 tt _1 r 

h 0, J 


dk' 


IV(k-k') 

1 


P k-k'> 2 ( 1 - cose kk' ) 6 ( V < Tc' ) >S 


( 2 ) 

with 

t, , = h 2 k ,2 /2m , £_ = h 2 k 2 /2m 

K r r 

and 

3 2 

K = 3 tt n 
F e 

Equation (1) represents the ensemble average of the resistivity 
calculated in Born approximation for elastic scattering from 
each configuration of the ions described by the density components 


P 


k-k' 


N 


-l 


-i(k-k'>R 


( 3 ) 


where ^R^ j- is the instantaneous set of ionic positions. The 
matrix element of the (self-consistently screened) aggregate 
electron ion scattering potential V(r) is defined for plane 
wave levels lk> by 

1 fsj 

fl < klvlk’ > = V(k-k’) = f dr e" 1 ^ - -' v(r) ( 4 ) 

n 
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If the scattering is sufficiently weak (Appendix A), Eqs. (1) 

3 

and (2) reduce, as originally shown by Ziman to 



where y = Ik-k'l^k-,, and v(y) is the electron-ion interaction 

J 1 rss 1 r 

scaled to its long wave length limit (2/3 ^p)* The quantity 
2 

a Q h/e may be viewed as the atomic unit of resistivity and has 
the convenient practical value of 21.7 uficm. S(k-k') is the 

iv 

liquid structure factor defined by 


S <3> - N ^ P q P -q ?> - 


N6 


( 6 ) 


In the Percus-Yevick model , (for hard spheres of diameter a), 
S(cp is a function of the packing fraction r\ given by 

(7) 

4,5 


r\ = T- n. CT 3 , n. = N/fl . 
1 6 ion ion 


For most classical fluids near their solidification points, 


r\ psi 0.45. 


We are dealing with point ions and the accuracy with which 
v(y) can be specified is limited only by the uncertainties in 
the dielectric function e (y) . In the neighborhood of y ~ 1 
(the regime dominating the integrand of (5)) e (y) is quite well 
known and we take the interpolation form suggested by J. 
Hubbard^ so that 

v(y) = -0.166 r / (y^ + 0.166 r F(y)) , (8) 

s s 


where 



- 6 - 


F(y) = f(y)/(l - 0.166 r f(y)(2y 2 + g)' 1 ) , 

b 

g = (1 + 0.0262 r ) _1 

s 

and 

f(y) = I + ‘ 

In practice, the replacement of F(y) in (8) by the Lindhard 
function f(y) leads to the same resistivity (to within 27 0 ), 
but the exchange and correlation corrections contained in F(y) 
are important in calculations of quantities involving (— - 1), 
such as the effective pair interaction between ions. 

3 Crr 

Since (r g a o kp) = we may rewrite (5) (using (8)) 
as 

p/(r 3 Z) = 38.4 \ dy y 3 S(y)(y 2 + 0.166 r F(y))" 2 
s J o s 

pftcm (9) 

The utility of this expression is that the right hand side 
is, for r < 1, a weak function of r and hence density. 

Figure 1 demonstrates this clearly. It is worth noting that 
the charge Z enters in the structure factor. ^ 

To obtain the resistivity as a function of temperature 
we require T(t\) at each density. This can be obtained from 

g 

a variational technique but the method is laborious and for 
the present purposes it is sufficient to use the approximate 
technique suggested by Ashcroft and Langreth^. We evaluate 
the pair interaction between point ions from 
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0(P) = (0.166 r g ) 


I 


xsxnpx 


dx 


o [x + 0.166 r F(x) ] 

s 


( 10 ) 


which gives the pair energy at separation r(r = p/2k ) in 

r 

units of £ (see Figure 2). If 0 . is the minimum value 
r mm 

of 0(p), then the melting temperature can be estimated 


from the relation 

V 


0(2^0) - 0 (min) = ( 3 / 2 ) k B T M /£ p 


provided 2k^a is evaluated at r\ = 0.45. It may be noted that 
this procedure gives in sodium to within 10%. To find 
dq/dT, we evaluate the slope of 0(p) 

d0 


dT " ( 2,l/P (dP ) T F ) p . 2 <lthZn) 1/3 


( 11 ) 


.5 ,2 


(where T F = 6 x 10 /r g K) and in this way obtain (see 
Fig. 3) and the values of T appropriate to q < 0.45. An 
alternative method for obtaining T^ exploits the Lindemann 
rule (see Appendix B), but the simpler approach outlined above 
is no less accurate and is, in fact, more fundamental. 

The results of our calculation for fully ionized H, He 
and C are found summarized in Figs. 4, 5 and 6. We choose as 
a vertical axis the quantity [resistivity x density] since 
as noted above this combination, near T^, is weakly density 
dependent. It should be emphasized that if our estimates of 
are incorrect, the form of the curves presented will remain 
substantially correct. We should also point out that at 
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densities for which the element carbon is likely to be fully 

pressure ionized, the hard- sphere approximation to the ion-ion 

9 

interaction may already depart substantially from reality. 
Moreover, kl^/£’ F ~ 0.05 and this implies a significant non- 
degeneracy. 

Fig. 7 shows a comparison of our results with those of 
2 

Hubbard and Lampe . The quantity compared is the conductive 
10 

opacity as tabluated in Ref. 2. Our results are seen to 
be systematically lower and the greatest difference occurs at 
low temperatures where the crude approximation for S(q) 
used in Ref. (2) is expected to be least accurate. We cannot, 
however, eliminate the possibility that the systematic dis- 
crepancy results from a disagreement in the temperature scale. 



III. Extension to Alloys 


The extension to binary alloys is straightforward in 


principle. The result equivalent to Eq. (5) can be written 

1 


(rV) 


= 38. 




3 . 

V dy 


o [y^ + 0.166 r F(y)]^ 

s 


{xS 22 (y) 


2x l (l-x) i S 12 (y) + (l-x)S n (y) } 

(lift- cm) . 


( 12 ) 


where x is the fractional number of ions of species 2, Z is 
the number of electrons per ion, and S 22 > are partial 

structure factors. ^ These structure factors not only depend 
on 

_ volume occupied by hard spheres 
^ total volume 


but also on 


a l 



where a^, 02 are the hard-sphere diameters of components 1 and 

2 . 


If a = 1, then Eq. (12) becomes identical to Eq. (9^ 

* 

except, of course, that Z is a function of x. In this special 

case, the results of Fig. 1 can be used to find the resistivity 
12 

of any alloy at the melting point. 

Eq. (10) shows that if 0(P q ) = 0 for the interaction between 
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ions of species 1 then 0(P q ) = 0 for the ions of species 2. 

i 

This suggests that a is near unity. However, the species with 

higher ionic charge is expected to have a "harder" core (for 

13 

a given value of r ). A detailed calculation suggests that 

s 

a - 0.75 for a hydrogen-helium mixture; that is, the helium 
hard- sphere diameter is one- third larger than the hydrogen 
hard- sphere diameter. In Fig. 8, we show that this deviation 
from a = 1 does not dramatically change the resistivity and 
accordingly a reasonable approximation sets all hard- sphere 
diameters equal. 

There is, however, no simple extension of our method for 
obtaining d-q/dT to the alloy problem. For Z > 2, the tempera- 
ture dependence of the resistivity is sufficiently weak that 
it may be ignored in a first approximation (for < T « T p ) . 
For a hydro gen- he lium alloy, a crude approximation simply 
interpolates between the temperature trends shown in Figs. 4 


and 5. 



IV. Summary and Application 


In the limited temperature and density range appropriate 

to Eq. (5) and the hard-sphere model, we find somewhat lower 

2 

resistivities than those previously obtained for fully ionized 
liquid metals. This is attributable to the use of a more 
accurate electron-ion interaction and a more appropriate 
structure factor. A disadvantage of the present method is 
the need independently to estimate the temperature scale. 

Systems for which the present calculations seem likely 
to apply include the interiors of the giant planets, in parti- 
cular Jupiter. Most recent models of the Jovian interior 
postulate a central region of dense fluid. Its composition 
is predominantly metallic hydrogen, but is augmented by a 
small amount of helium (about 107 o by number^*^) . it is 
conceivable that the helium may not be completely ionized 
and if not, the electron-helium interaction may be more appro- 

■I £ 

priate to that expected of neutral helium atoms. We find that 
although it is possible for the resistivity to be enhanced 
if the helium remains un-ionized, this enhancement is mainly 
a consequence of the small increase in the value of r g rather 
than any substantial change in the scattering cross section 
from that expected for fully ionized atoms. 

If we choose the central temperature^ of Jupiter to be 
about 16,000°K, then we find that the resistivity of the fluid 
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is expected to range from 4 pft-cm at the center of Jupiter 
to about 8 (jfi-cm at the boundary between metallic and molecular 
hydrogen. A conductivity characteristic of the deep interior 
of Jupiter is therefore 
a ~ 2 x 10^ esu 

18 

a result somewhat larger than most previous estimates. 

Jupiter is observed to have a strong magnetic field and 
in seeking internalmechanisms for its origin it is first of 
interest to decide whether the field could be primordial. 

If it were, then the quantity of central importance is the 
decay time T given in seconds by 
T ~ 4 tt a (L/ c) ^ 

where c is the velocity of light and L is a typical planetary 

9 

dimension, which we take here as 5 x 10 cm. The result 

9 

T ~ 2 x 10 years 

may seem to hinge not too seriously on the choice of L. Even 

if the value chosen is viewed as unreasonably large, the 

result for T remains such that the possibility of primordial 

origin is difficult to discount. In complete contrast to this, 

it is interesting to record that the high value of ct is likely 

19 

to be favorable for a dynamo mechanism underlying the generation 
of the magnetic field. 

Finally, a straightforward application of the Wiedemann-Franz 
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relation yields thermal conductivities for the interior of 

g 

Jupiter ranging from (in erg/cm-sec-K) 9 x 10 at the center 
8 

to 1 x 10 at the metallic boundary. Now the observed internal 

20 

heat flux is very high , but it is apparent that even con- 
ductivities of this magnitude are insufficient to maintain 

the measured flux unless we assume a much larger central 
14 15 

temperature * In a situation such as this, the system 
is unstable against convection and the planet would rapidly 
cool. It would seem to follow that all but a small core of 
Jupiter must be convective. The size of this convective 
region is an open question. 
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Appendix A : Validity of the Born Approximation 

An elementary criterion for the validity of the Born 

2 

approximation is that 
„ 2 2 h 

aa_ < (J!_y . 

e p ~ V2me F > 

Here, the left-hand side is roughly the distance from the ion 
within which the interaction energy exceeds the Fermi energy. 
The right-hand side is of the order of the electron wavelength. 
It follows that 
2 2 

27 e 2 

£-, > o — = 4Z Rydbergs 

(-*?) 

me 


whence r < — . 
s - Z 

An alternative criterion is 


Born 


« 1 


4tt a * 


2 . 


where 4na is the "geometric" cross section. For a single ion 

2k F _ _2 


Born 


= 7T2 J *( 5 v(k > ] * 


2rtkp v o “ h 


where 


V(k) = 1* sinkr V(r) rdr. 

o 

We calculate ^g orn approximately using Thomas-Fermi screening. 


i.e 


= Ze -q s r 


V(r) = =j-e 
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so 


2 2 
„ 4TTZe _ 4rrZe 

v(k) = TTTT ‘ T 


k +q x +0.166 r (2k )‘ 


where 


x = 


2k F ' 


Thus, 


„ - 1 1 

r4nZ \ 

z . 1 1 

r xdx 

°Born _ . 2 

2uk F 

^ a o ^ 

(2k F ) 2 J 

o (x 2 +0.166 r ) 2 
s 


But, 


a ~ — = 0.64 r^a 
— q so 


and thus it follows that 

2 2 
0.27 r Z z 
Born s 


4tra 2 - < 1+0 - 166 r s> ' 


Finally, 

« 1 implies r < \ (as before). 

. Z. S f*> L 

4ua 

However, the Born cross section per ion in the condensed 
state is clearly different from that of a single isolated ion. 
We can calculate the "apparent n cross section, per ion, in 
the liquid by using the identity 
no a v F T e 1 

where t is the "collision time" for an electron and n is the 

ion number density. 

m 


2 

n e t 
el 


Since P = 
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we have, from Eq. (5) 


o = 4]L 2 - Z - J y 3 v 2 (y) S (y)dy 
k„ o 




whence 


a 


4rra' 




where P is calculated from the first Born approximation. 

(Note that this formula is valid for any simple liquid metal.) 
For hydrogen at r = 1.6, T = T^, 


4rra 


2 - 


0.06, 


and for helium at r = 1.2, T - 1^, 


a 


4na 


2 - 


0.25. 


This suggests (but does not prove) that the Born approxi- 
mation may be much better satisfied in the condensed state 
than for a single ion. Thus, our criterion r < ~ may be too 
stringent. It is clear and expected, however, that the Born 

approximation is increasingly well satisfied as r becomes 

s 


smaller. 




where y is the mean square amplitude of the ions just below 
the melting point and is found, almost universally, to be about 
1/16. M is the ion mass, R q the interatomic spacing, 
a phonon frequency of wavevector k and polarization X and 
is the Bose-Einstein occupation factor. 

22 

For the high-density systems considered, Abrikosov 
has shown that it is important to distinguish between the 
longitudinal and transverse modes, since the former are primarily 
determined by the bulk compressibility of the electron gas 
whilst the latter are primarily determined by the Coulomb 
forces between ions. 


We make a Debye approximation but allow for the longi- 
tudinal and transverse "Debye" temperatures to be different. 

23 

Using the method outlined by Trubitsyn we obtain (in K) 



- 17 - 



-18- 


(the correlation energy of the electron gas is small and can 


be ignored). Eq. (Bl) can then be written 

s./t 


MS. t o 


l 

Vt 

MS 2 


X 1 
e - 1 


] 


0 /t 




xdx 


x i 

e - 1 


] ~ 0.47 


(B2) 


where S^, S t are the appropriate sound velocities. We anti- 
cipate T^ < 0^, 0^ and so approximate 0^/T, 6^,/T by <» in the 

integrals. It is easy to show that this is valid provided 

2 

vj; e for 0=0 and 0 , 

T 6 i t 

which is satisfied reasonably well for the cases studied. 

Eq. (B2) can be written in numerical form, for low temperatures, 


as 


0.22 


A i z 5 /6 (2^1 


3.66 


2/3 > [ 1+ ¥(f) 2 ] 

7.17 Z 2/3 \* . 2 L 3 V0 X ' J 


) 


+ ot [ 1 + ^)'> 0 - 47 


2 2 
2tt / T ‘ Z 


and is solved to obtain T^ (note that for r g « 1, only the 
transverse modes are important in determining T^) . The results, 
shown in Fig. 9, give melting temperatures which differ by as 
much as a factor of two from those in Fig. 3. Similar results 
have been obtained by Pollock and Hansen. The problem 
with Lindemann's rule is that an error in y (= 1/16 in the 
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above calculation) propagates alarmingly through to the final 
calculation of T^, in the case T^ < 9 t » 0^. Typically, a 
107« error in y will give a 50% error in T^. Moreover, our 
estimates of 0 , 0 are only approximate. (Our formula for 

u Xj 

0 is, however, in excellent agreement with the 9_ calculated 

25 

by Neece, Rogers and Hoover. ) Note that at sufficiently 
high densities, the zero point motion alone will cause the 
lattice to melt. Lindemann's rule gives an estimate of the 
value of r® at which 1^ -* 0. Since density varies as (r^) > 

the density at which -♦ 0 cannot be calculated to better 
than an order of magnitude using Lindemann's rule. As 
Abrikosov‘S observes, only hydrogen and helium will melt at 
absolute zero and sufficiently high densities. This is because 
the densities required for heavier elements are such that the 
sizes of the nuclei become important. 


The pressure at which T^ - 0 may be incorrect by 
almost two orders of magnitude. 
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Resistivity of fully-ionized liquids at r\ = 0.45. 
Effective ion-ion interaction energy in units of 



Estimated melting temperatures. 

Resistivity of hydrogen. 

Resistivity of helium. 

Resistivity of carbon. 

Conductive opacity of hydrogen at two densities. 

A comparison of our results with those of Hubbard and 
Lampe (Ref. 2). 

Resistivity of an H-He alloy at r =1.0 and q = 0.45. 

s 

The effect of different hard sphere diameters is shown. 
The melting temperatures of metallic hydrogen and 
helium according to Lindemann's rule. 
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